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ABSTRACT 

Recent observations by PAMELA, Fermi-LAT, and AMS-02 have conclusively indicated a rise in the 
cosmic-ray positron fraction above 10 GeV, a feature which is impossible to mimic under the paradigm 
of secondary positron production with self-consistent Galactic cosmic-ray propagation models. A 
leading explanation for the positron fraction rise is an additional source of electron-positron pairs, for 
example one or more mature, energetic, and relatively nearby pulsars. We point out that any one of 
two well-known nearby pulsars, Geminga and Monogem, can satisfactorily provide enough positrons 
to reproduce AMS-02 observations. A smoking-gun signature of this scenario is an anisotropy in the 
arrival direction of the cosmic-ray electrons and positrons, which may be detectable by existing, or 
future, telescopes. The predicted anisotropy level is, at present, consistent with limits from Fermi- 
LAT and AMS-02. We argue that the large collecting area of Atmospheric Cherenkov Telescopes 
(ACTs) makes them optimal tools for detecting such an anisotropy. Specifically, we show that much 
of the proton and 7-ray background which affects measurements of the cosmic-ray electron-positron 
spectrum with ACTs may be controlled in the search for anisotropics. We conclude that observations 
using archival ACT data could already constrain or substantiate the pulsar origin of the positron 
anomaly, while upcoming instruments (such as the Cherenkov Telescope Array) will provide strong 
constraints on the source of the rising positron fraction. 



1. INTRODUCTION 

Local observations of cosmic-ray fluxes are among the 
most powerful probes of high-energy astrophysical pro- 
cesses in our Galaxy. These measurements can probe 
both the nature of particle production in energetic as- 
trophysical systems, as well as the physics of particle 
diffusion through the Galaxy. Recently, results from the 
Advanced T hin Ionization Cal orimeter (ATIC) balloon 
experiment ( Chang et al.|200^ ), as well as from the Pay- 
load for Antimatter Matter Explor ation and Light-nucle i 
Astrophysics (PAMELA) satellite ( [Adriani et aLl|2009[ ), 
indicated possible anomalies in the expected flux of cos- 
mic ray leptons. Specifically, ATIC observed a bump in 
the cosmic-ray electron-plus-positron (e^) spectrum at 
an energy of approximately 800 GeV, the intensity of 
which has since been cast into doubt by observations 
with the Hig h Energy Stereosco pic System (H.E.S.S.) 
telescope ( [Aharonian et al. 20091) and the Ferrni Larg e 
Area Telescope (l-'erml-LA'i ji Ackermann et al. 2010a). 
PAMELA, on the other hand, using a built-m niagnetic 
field that allows it to differentiate the charge of cosmic- 
ray species, observed a rise in the positron fraction of 
the cosmic-ray e* spectrum above 10 GeV. This result 
has been confirmed by subsequent observations using 
the Fermi-LAT, cleverly employing the geomagnetic field 
structure and the Earth shadow to disting uish cosmic-ray 
particle charges ( Ackermann et aL]|2012 ). 

Very recently,^ tEe Alpha Magnetic Spectrometer 
(AMS -02) released its initial findings (Aguilar et al. 
20131, which confirm the previously observed rising 
positron fraction with significantly better statistical and 
systematic errors than previous experiments. Addition- 
ally, AMS-02 extended the energy range of the positron 
fraction up to an energy of 350 GeV, finding that while 



the positron fraction continued to rise at these ener- 
gies, the slope of the positron fraction appears to flat- 
ten above ~150 GeV. This suggests that the positron 
fraction may either approach a stable asymptotic value, 
or that it might even decrease at higher energies. How- 
ever, no evidence of a sharp cutoff in the positron frac- 
tion was observed. Compared with PAMELA measure- 
ments, AMS-02 results also find an anisotropy which is 
approximately 20% smaller at high energy. Since the sec- 
ondary positron signal in this range is extremely small, 
this results in an approximately 20% reduction in the 
primary positron flux necessary to explain the AMS-02 
positron ratio in light of Fermi-LAT constraints on the 
total electron plus positron spectrum compared to fits to 
the PAMELA measurements. In addition, the slightly 
softer spectrum inferred from the AMS-02 data com- 
pared to the PAMELA measurement also points to a 
smaller normalization of the needed additional primary 
component. 

Myriad models have been formulated in order to ex- 
plain the rise in the positron fraction, nearly all of which 
require a new injectio n spectrum of primary positrons ( |di] 
Bernardo et al.|2011 ). These models can be broken down 
into two major categories. The majority (too numer- 
ous to list here) concern the annihilation of dark mat- 
ter particles into, primarily if not exclu s ively, l eptonic 
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final states (e.g. [A rkani-Hamcd et al. (2009); Cholis 



et al.| ( |2009| ); |CirellTet al.| ( ,20U 9)). RelevanTHafk matter 
models almost uniformly require a massive dark mat- 
ter particle (M^^ > 0.5 TeV, possibly even larger af- 
ter the AMS-02 data extended the positron fraction to 
350 GeV) with an annihilation cross-section that signifi- 
cantly exceeds the thermal annihilation cross-section im- 
plied by the "WIMP miracle". Since leptophilic dark 



matter annihilation should additionally create a bright 
7-ray signal (for example from inverse Compton up- 
scattering of background radiation fields, or from final 
state bremsstrahlung, or from the hadronic decays of 
T leptons), many of these models have been ruled out 
by Fermi-LAT 7-ray observations of, e.g., the Galactic 
center, dwar f spheroidal galaxies or nearby galaxy clus- 



ters (see e.g. Meade et al. 



2010 Hooper & Linden 2011 



Hooper et al. 2012[ |Ackermann et al. 2011; Acke rmann 



et a .| 2010[ for early studies). Other constraints arise 
trom distortions to the CMB spectrum or from radio 
observations. These analyses, will need to be updated 
in view of the new AMS-02 data, which surely restrict 
the class of dark matter models which can explain the 
positron fraction solely based on the observed positron 
fraction spectral shape. However, this task lies beyond 
the scope of the present study. 

A second class of models relies o n mature (i.e. in a n 
age range between lO"* and 10^ yrs ( Grasso et al.||2009 )) 



pulsars, which are known sources of e pairs, in order to 
produce the rising positron fraction. Models of particle 
creation in pulsars' magnetospheres indicate that these 
objects produce a significa nt population of energetic lep- 
tons with a hard spectrum ([Hooper et al.|2009 Delahaye 
et al."201 0[|Profumo||2012||Maiysiiev et al.||!^0o01|iJarger 
et al. 2009| [Kistler &: Yuksell feoOQp . The pulsar con- 
tribution to the positron traction could stem from the 
cor abined fluxes of a G alactic population of pulsars (see 
e.g. Barger et al. 2009|), or may be dominated by a sing le, 
nearby source (eg. |Profumo||2012[ [Yliksel et aL]|2009[ ). 

Pinpointing the origin ot the rising positron traction re- 
quires a method of differentiating pulsar and dark mat- 
ter models for this signal. While each class of models 
may produce slightly different predictions for high en- 
ergy diffuse radiation (e.g. from inverse Compton emis- 
sion), our poor knowledge of the diffuse 7-ray sky has 
precluded this avenue from conclusively distinguishing 
between models. Similarly, accurate measurements of 
the total e* spectrum can hardly produce a smoking 
gun favoring either model: in the dark matter case, a 
spectral cutoff in the positron fraction emerges at the 
kinematic edge corresponding to the dark matter par- 
ticle mass; in the case of pulsars, a cutoff is predicted 
from energy losses and diffusion and/or from the intrin- 
sic source injection spectrum. An observational test that 
could conclusively rule out a dark matter interpretation 
does, however, exist; the detection of an anisotropy in 
the arrival direction of cosmic-ray leptons. 

Models using dark matter to explain the cosmic-ray 
excess are unlikely to produce significant anisotropics in 
the arrival direction of the e^ observed at the solar posi- 
tion. A dark matter clump that produced a large enough 
anisotropy would also produce a bright gamma-ray emis- 
sion that would be easily and disti nctively detectable by 
the Fermi-LAT (|Brun et al.||2009|). The anisotropy in- 



duced by nearby pulsars m ay instead be significant and 
potentially observable (Hoo per et al.| 2009), although the 
anisotropy may be highly suppressed if numerous pul- 
sars located randomly in the s ky contributed si gnificantly 
to the overall positron flux. Profumo (20121, however, 
pointed out that the expected anisotrojaies from selected 
nearby pulsars (e.g. Geminga, Monogem, Vela, CTAl, 
B0355-I-54) may be detectable at a 2a level with 5 years 



of Fermi observations. A search using one year of F ermi 
LAT data did not un cover any signs of anisotropy ( Ack 
ermann et al. [[20^1 Ob ), consistent with expectations from 
pulsar models. Intact, the current constraints lie approx- 
imately an order of magnitude above pulsar projections. 

Here, we ague that Atmospheric Cherenkov Telescopes 
(ACTs) offer a very promising tool in searches for cosmic- 
ray e* anisotropies associated with a local source. Due to 
their large effective area, ACTs are capable of detecting 
much smaller anisotropies than space-based telescopes. 
However, ACTs also face significant hurdles due to the 
systematics associated with differentiating e* showers 
from hadronic showers created by high energy protons, 
which dominate the e^ flux by several orders of magni- 
tude: Uncertainties in hadronic shower rejection present 
a systematic error of approximately 30% in the cosmic- 
ray e* flux observed by H.E.S.S., an uncertainty which 
is several times larger than those from the Fermi-LAT 
cosmic-ray e* measurement at similar energies. 

In this paper^ we first point out that a simple model 
employing a single nearby pulsar in addition to the best- 
fit diffuse cosmic-ray background produces an extraor- 
dinarily good fit to the positron fraction observed by 
AMS-02 and is consistent with the combined electron 
plus positron flux as measured by the Fermi-LAT and by 
H.E.S.S. The AMS-02 data are perfectly compatible with 
estimates of the e* flux from nearby mature pulsars such 
as Geminga and Monogem, although other pulsars may 
also contribute. We note that these pulsars are expected 
to produce a significant anisotropy in the cosmic-ray e* 
spectrum, and evaluate the associated signal. While 
this anisotropy falls nearly an order of magnitude be- 
low the current constraints from both AMS-02 and the 
Fermi-LAT, we point out that the large effective area of 
ACTs may produce significantly more robust constraints. 
Specifically, we show that as long as the hadronic back- 
ground observed by these telescopes is highly isotropic, 
it can be treated as a statistical, rather than a system- 
atic uncertainty. Already existing ACT observations are 
likely to significantly cut into the scenarios where, for ex- 
ample, the Monogem pulsar is responsible for the rising 
positron fraction. Finally, we note that the upcoming 
Cherenkov Telescope Array (CTA) will greatly enhance 
the potential for searches for cosmic-ray e^ anisotropies, 
due to both its improved effective area and hadronic re- 
jection efficiency. 

The outline of this paper is as follows: In Section [2J 
we discuss the contributions of young, known pulsars to 
the total electron plus positron flux observed at the solar 
position. In Section [3] we compute the anisotropy result- 
ing from the Geminga and Monogem pulsars which we 
find to provide the best fits to the observed cosmic-ray 
signal, and we discuss the possibility for detecting this 
anisotropy with ACTs. In Section HI we present our re- 
sults on anisotropy searches with ACTs, and discuss pos- 
sible systematics resulting from Galactic diffuse emission 
in Section [5] Finally, in Section [7] we discuss the impli- 
cations of our results and conclude. 

2. THE CONTRIBUTION OF NEARBY PULSARS TO THE 
DIFFUSE ELECTRON-POSITRON FLUX 

In order to calculate the background cosmic-ray elec- 
tron and positron flux attributable to diffuse Galactic 
processes, we employ the Galprop cosmic-ray propaga- 



tion code (Strong & Moskalenko 1998 Strong et al. 
2009 1 . Galprop selt-consistently calculates the injection 
and diffusion of cosmic-ray electrons, positrons, protons, 
and nuclei throughout the interstellar medium of the 
Milky Way galaxy, taking into account effects such as 
convection and re-acceleration, as well as energy loss pro- 
cesses such as synchrotron and bremsstrahlung emission. 
The Galprop code has been specifically tailored to ac- 
curately reproduce the ratios of radioactive cosmic-ray 
primary and secondary nuclei (e.g. the Boron to Carbon 
ratio) which currently stand as the best handles on the 
parameters of cosmic-ray propagation in the Milky Way. 
Specifical ly, we employ here the best-fitting Bayesian 
model of Trotta et al. (20111, which finds a posterior 
mean for the diffusion constant of 8.32x10^^ cm^s~^, 
normalized at 4 GeV, with an energy index S = 0.3, 



nearly equivalent to that first proposed by [Kolmogorov 
(1941). We employ all parameters from this model in 



the following analysis, but linearly vary the total injec- 
tion normalization of cosmic-rays in order to create a best 
fitting model when various pulsar profiles are added to 
account for additional sources of electrons and positron^ 
We calculate the contribution to high-energy cosmic- 
ray electrons and positrons from pulsars in the co ntext 
of the point-li ke and burst- like approximation of |Shen| 
( 1970 1 (see also Panov ( 2013 1 for a recent review). We as- 
sume that pulsars mject m tlie interstellar medium (ISM) 
as many electrons as positrons, with a source function 



Qif,t,E)^Q{E))S^r)6it)), 



(1) 



with an exponentially suppressed power-law spectrum 
parameterized by 



QiE) = QoE-'e^p{E/E,^t). 



(2) 



The normalization Qq of the source function is then as- 
sociated with the total energy output of the pulsar via 
the simple relation 



Q{E)EAE = TjWo, 



(3) 



with Wo indicating the total pulsar energy output, and 
77 < 1 the relative efficiency with which such energy is 
converted into the energy associated with e* pairs in- 
jected in the ISM. 

The source function of Eq. (fTl) is substituted in the rel- 
evant differential equation describing cosmic-ray trans- 
port in the limit where the energy loss terms and the dif- 
fusion constant are both independent of positio n, yield- 



ing a simpl e, exact analytic solution (see e.g. Atoyan 
et al. 1995). Among the remarkable features of this 



asymptotic solution is that there exists a cut-off in energy 
set by the e inverse Compton and Synchrotron energy 
losses and by the time T of the burst-like injection. If 
the energy loss term is cast as 



d^ 

dT 



= -ba E\ 



(4) 



the cutoff corresponds to E\oss — ^/{boT). Here, we em- 
ploy bo = 1.4 X 10^16 GeV"^s-i. This implies Ei^ss ^ 
750 GeV for T = 300 kyr, and Eioss ^ 226 GeV for T = 1 



^ In practice, wo adopt a simple rescaling of the Galprop diflfuso 
e^ flux by a factor 0.8. 



Myr. On the other hand, there exists a second exponen- 
tial suppression, due to e^ diffusion during the relevant 
energy-loss time-scale, of the form 



exp 



{d/riis{E)y 



(5) 



where d is the distance to the (point-like) source and 
TdieiE) is a characteristic, energy-dependent diffusion 
length. For example, at energies around 1 TeV the typ- 
ical diffusion length falls well within 1 kpc. These two 
constraints effectively limit the class of pulsars that could 
add to the positrons feeding the anomaly observed in 
the rising positron fraction. In particular, the age of the 
candidate pulsar must fall between the time-scale for the 
confinement of e pairs within the nebula (~ 10* yr), and 
the time-scale corresponding to energies large enough to 
explain the highest energy where the positron anomaly is 
observed {^ 1 Myr for E ^ 350 GeV ). Pulsars in this ag e 
range are usually dubbed "mature" ( Grasso et al.|2009 ). 
While we do not attempt here to exhaustively pro- 
vide a census of the known pulsars that could contribute 
to the positron fraction anomaly, this exercise has been 
carried out in the p ast for the PAMELA data (e.g. |Pro 
fumo (2012), Grasso et al] (|2009[) , [Malyshev et al. (2009 
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Hooper et al., (|2009j) among others). Instead, we provide 
two examples of well-known pulsars in the right distance 
and age range to potentially contribute very significantly 
to (or as we shall see, explain com pletely!) the recent 
high-statistics AMS-02 data ( jAguilar et al.||2013| ). 



We focus here on two of the most luminous nearby pul- 
sars in the correct age and distance range: the Geminga 
(J0633-t-1746) pulsar and the Monogem (B0656-I-14) pul- 
saiFl We assume that the burst-like injection happened 
for Doth pulsars at look-back times comparable to the 
age of the pulsar (i.e. we neglect the trapping time of 
e* in th e nebula compared t o the age of the pul- 
see e.g. 



ip; 

3' 



Grasso et al.| (2009)). We set Geminga's 



10*^ yrs, and Monogem's to 1.11 x 10 yrs. 



the 
sar 

age to 3.42 x 

and the distanc es to each pulsa r to 0.15 and 0.29 kpc, 
respectively (see Profumo]|2012 and references therein). 
While the age and distance for the pulsars is set by ob- 
servation, less is known about the spectral shape and en- 
ergy output. The exponential cutoff -Ecut is entirely im- 
material for Geminga, as the corresponding i^ioss ^ ^-cut, 
while for the purpose of reproducing the AMS-02 results 
Ecut is also irrelevant for Monogem, unless it falls in 
the sub- TeV domain, which is theoretically implausible 
(Malyshev et al. 2009). For definiteness, we set for both 



pulsars i/'cut = 2 I'eV (changing this parameter clearly 
affects the details of the spectrum of e^ a.t E > Scut, 
but reproducing such spectrum in detail is not within 
the scopes of this study, as it is largely affected by the 
diffuse Galactic background, not accurately modeled by 
Galprop at these energies). 

Direct observations of the pulsar source spectral index 
7 are available for a handful of pulsar wind nebulae (see 
e.g. [Green 2009) , and point to 7 < 2, althoug h in some 
cases 7 ~ 2 is also observed (Aharonian et al. 2006). 



We find that, for our choices of diffusion parameters, the 
AMS-02 data are best fit by 1.8 < 7 < 2, and we em- 

^ Other radio (or 7-ray only) pulsars are also capable of con- 
tributing significantly to the local positron f lux. A more thorough 
discussion of such candidates can be found inlGendelev et al.|l|2010| 
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Fig. 1. — Left: The positron fraction from a combination of the Galprop model for the diffuse e^ Galactic background (green dotted), 
along with contributions from the Geminga (black) and Monogem (red) pulsars, compared with data from PAMELA (green circles), Fermi- 
LAT (orange triangles) and AMS-02 (blue squares). Right: The flux of cosmic-ray electrons and positrons from a combination of the same 
Galprop model (green dotted), with contributions from the Geminga (black dashed) and Monogem (red dashed) pulsars. These create 
a total cosmic-ray lepton spectrum (black and red solid respectively), which can be compared with data from the Fermi-LAT (orange 
squares) and H.E.S.S. (pink diamond) observations, (right). Note that the diffuse background from Galprop was not tuned to reproduced 
the H.E.S.S. data, and we do not attempt to fit those data above 1 TeV. 



ploy 7 = 1.9 for Geminga and 7 — 1.95 for Monogem. 
The resulting normalizations required to provide a fit 
to the AMS-02 data with a single pulsar correspond to 
r]Wo = 2 X lO''^ erg for Geminga and to 77 Wq — 8.6 x 10^^ 
erg for Monogem. Within the context of our diffusion 
model, we note that these values act as upper limits on 
the total lepton flux from each pulsar for any scenario 
which is compatible with the AMS-02 data, since these 
values must decrease if additional sources are considered. 
The total energy outputs we find depend quite sensi- 
tively on the assumptions made for the spectral slope, 
but are generically compatible with the total energy out- 
put expected from a mature pulsa r, which ranges within . 



5 X 10^^ < Wn/eiK < 5 x 10^°, (Delahaye et al. 



Malyshevet al. 20091. 



i 



2010 



Employing a combination of the Galprop Galactic e 
diffuse background model, rescaled by a factor 0.8 to 
account for the additional sources, and the calculated 
flux from each candidate pulsar, in Figure [T] we show the 
positron fraction (left) and the combined flux of electrons 
and positrons (right) observed at the solar position for 
models in which the Geminga pulsar dominates the pro- 
duction of nearby positrons (black), and a model where 
the Monogem pulsar dominates cosmic-ray positron pro- 
duction (red). In each case, we find an extremely good 
match between our results and AMS-02 observations. 

3. DETECTION OF A COSMIC-RAY 
ELECTRON/POSITRON ANISOTROPY WITH ACTS 

In the context of diffusive propagation, we estimate the 
expected anisotropy from a source at a d istance d that 
injected e^ at a time T (e.g. Grasso et al.||2009| ) with 



A = 



3 d {l-6)E/Ei,,, NpsAE) 



(6) 



2cTi_(i_ij/i5,_)i-^Artot(i?)' 
and A'tnt the pulsar and total e* spectra. The 



with A''psr aiiu ivtot 

dipolar anisotropy A is defined as 



A = 



Nf + Nb 



(7) 



where Nf and Nj, are the total number of e* ob- 
served during a selected ensemble of observations point- 



ing within the sky hemisphere centered on the pulsar 
(Nf) and during a second ensemble of observations with 
the same collective effective exposure as the first ensem- 
ble, pointing within the opposite hemisphere (iVh). 

It is worth noting that this calculation of the 
anisotropy from a single pulsar is overly simplistic, as ig- 
nores several possible complicating effects. For instance, 
the corresponding anisotropy might be washed out by ef- 
fects such as a local magnetic field bubble, the pulsar's 
proper motion during the age of e* injection, or signif- 
icant deviations from the simple diffusive propagation 
setup employed to theoretically estimate the anisotropy 
( Profumo||2012 ). On the other hand, anisotropics in the 
charged cosmic-ray spectrum can also be induced via dif- 
fusion in the interstellar me dium, for instance by local - 



magnetic field anisot ropics ( |Drury fc Aharonian 2008 
Giacinti fc Sigl|[2012 ). While this may produce a spun 
fon 01 



ous detection ot an electron/positron anisotropy not due 
to a nearby primary source, the two effects may be in 
principle disentangled in the following ways. First, any 
anisotropy induced by anisotropic diffusion should affect 
protons and electrons similarly, leading to a strong cor- 
relation between observed anisotropics for both species. 
In the case of a nearby e~^e~ source, which would not 
produce many protons due to the strong constraints 
on primary anti-proton production, the morphology of 
the anisotropy would not be seen in relativistic pro- 
tons. Second, any anisotropy stemming from particle 
diffusion is likely to have an anisotropy which depends 
on the scale of the magnetic field inhomogenities, while 
the electron anisotropy from a nearby source will have 
an energy dependent anisotropy which scales with the 
positron fraction due to that source. In particular, the 
anisotropy should disappear above any cutoff energy the 
primary positron source would possess. Lastly, inho- 
mogenitie s in diffusion parameters are likely to appear as 
hotspots ( |Drury fc Aharonian|2008{ [Giacinti k Sigl|20l2 1 
or streams! Kistler et al.|2012| in tiie data, an anisotropic 
signature from which is distinct from the dipole domi- 
nated term stemming from nearby sources. 

We now turn to the question of how to search for an 
anisotropy in the cosmic-ray e^ flux with ACTs. The 



most significant uncertainty in the determination of tlie 
cosmic-ray e^ spectrum with. ACTs is the efficiency of 
cosmic-ray proton rejection. This is the dominant sys- 
tematic error because the fiux of cosmic-ray hadrons 
dominates the lepton fiux by several orders of magni- 
tude. While observations of 7-ray point sources are able 
to employ the isotropy of the cosmic-ray signal in order 
to control this background, measurements of the cosmic- 
ray e flux must instead determine the hadronic or elec- 
tromagnetic nature of each individual observed shower. 
To this end, the H.E.S. S. collaboration has adopted a 
random for est approach (Breiman & Cutler 2004 Bock 
et al. 2004) intended to convert information about the 
observed shower into a parameter ^ which describes the 
extent to which the shower is electron-like. The param- 
eter C is determined in the range of - 1, with larger 
numbers indicating a better fit to Monte Carlo models 
of electron showers. While the C parameter is highly 
energy-dependent, in many situations its discriminating 
power is significant enough to produce an electron popu- 
lation which dominates the hadronic background at high 
C values. We note that even for moderate values of (, the 
contribution from heavier nuclei is entirely negligible. 

While a proper selection of C is important so that the 
cosmic-ray e* population produces a reasonable portion 
of the total cosmic-ray signal, searches for anisotropy 
are significantly less affected by errors in hadron rejec- 
tion compared to measuremen ts of the total e* spectrum 



([Aharonian ct al.'"2008', |2009| ). Assuming that both the 
misidentihed cosmic-ray proton and background cosmic- 
ray e^ fluxes are isotropic, the fraction of the background 
which stems from each is irrelevant in searches for e^ 
anisotropics. Instead, the measurable quantity is the 
fraction of the total cosmic-ray flux (with ( > 0.9) which 
stems from an anisotropic candidate pulsar. We can cal- 
culate the total contribution to the detected cosmic-ray 
flux with an ACT as: 



Ntot = (Npsr + N^) + iNe,,so + Np) 



PJt 



(8) 



where Npsr is the number of cosmic-ray leptons produced 
by the pulsar, Nj is the number of 7-rays observed as 
electromagnetic showers in the instrument, Ng^iso is the 
number of e from sources which are highly isotropized 
before diffusing to the solar position, and Np is the num- 
ber of protons misidentified as electromagnetic showers. 
We note that sources in the first parenthesis are mea- 
surably anisotropic, while the second set of sources are 
assumed to be perfectly isotropic in what follows. In 
our formal calculation, we will ignore the parameter Nj, 
noting that the 7-ray flux is highly subdominant to the 
e and proton fluxes. However, it is important to note 
that some 7-ray components, such as the Galactic diffuse 
emission, are a potentially important source of system- 
atic errors, because their anisotropy is of order unity. 
We discuss this systematic effect and the impact of our 
assumption at length in Section [51 

If N^ ~ 0, the dipole anisotropy can be written in 
terms of the pulsar contribution to the total cosmic-ray 
flux as: 

Nf + Nb Npsr J + Npsr,b + '2{Ne,rso + Np)' 

where, consistently with the notation introduced above 



in Eq. (|7|), Npsrj and Npsr,b are the total number of 
e* from the candidate pulsar from the two ensembles of 
observations when the telescope is pointing at locations 
in the hemisphere oriented directly towards, or directly 
away from the pulsar (note that the factor 2 in the de- 
nominator stems from our assumption of identical effec- 
tive exposures for the two ensembles of observations). 
An anisotropy measurement is then significant at the 2cr 
level when: 



A > 2 



N 



(10) 



avg 



where Navg = {Nf + Nb)/2 is the average number of 
showers observed during each ensemble of observations. 
Using this model, we can now estimate the anisotropy 
level, as a function of energy, which can be detected with 
an ACT, and compare it with the current observational 
constraints from Fcrmi-LAT and AMS-02. 

The key input to our estimate of the sensitivity to an 
e anisotropy is the determination of Navg which, in 
turn, depends on the e^ flux measured by the ACTs. 
We start with predictions on the limits which can be set 
using archival data from the H.E.S.S. telescope. We note 
that other ACTs, such as VERITAS ( [Holder et alpO OS ) , 
can in principle make very similar measurements. How- 
ever, details of VERITAS measurements of the electron 
spectrum are not publicly available, and we therefore em- 
ploy H.E.S.S. for our case study. 

H.E.S.S. has presented two different results on the 
cosmic-ray e^ spectrum. The first result conce rns the 



high-energy spec trum, in the range of 1-4 TeV (Aharo 



nian et al. 2008). In this regime, H.E.S.S. reports an 
effective area of approximately 5x10* m^. The study 
employed 239h of live-time, producing an effective expo- 
sure of 8.5x10^ m^ sr s. In this study, H.E.S.S. found a 
best-fitting e^ spectrum 



dA^/d£; = k{E/l TeV)" 



(11) 



3.9±0.1 (stat). Given the H.E.S.S. effective exposure, 
this implies the detection of 3370 e*. Indeed, using the 
Dexter packag^ we find that H.E.S.S. observed 3100 e^ 
with ( > 0.6. Of these events, 2600 were found to have 
C > 0.9, giving this cut a relative acceptance of 0.77. 
Interestingly, while H.E.S.S. reports the detection of 7610 
protons with C > 0.6, only 2470 of these protons have 
C > 0.9, yielding a sample which is, therefore, electron 
dominated. 

Additiona lly, H.E.S.S. investigate d the low energy e"^ 
spectrum in Aharonian et al. (2009), in ord er to investi- 
gate the spectral bump observed by ATIC ( Chang et al. 
20081). In this study, H.E.S.S. employed only 77 hours 



of livetime, with a calculated effective area of 3x 10 m , 
for a total effective exposure of 2.2x10^ m^sr s. H.E.S.S. 
obtained a best-fitting broken power law spectrum: 

dN/dE = k{E/Eb)-'^'{l + (S/^fc)i/")-(r^~ri)°, (12) 

with a best fit values Eb^O.9 ± 0.1 TeV, 
k = (1.5 ± 0.1) X 10-* TeV-^m-^sr-is-i, ri=3.0±0.1 
and r2=4.1±0.3. This indicates that H.E.S.S. observed 
approximately 7660 e^ with energies between 340-700 

^ http://dcxter.sourceforge.net/ 



6 



GeV, the range in which they provide the observed 
distribution over (.. Unfortunately, in this low energy 
regime, the C parameter is less efficient at differentiating 
the cosmic-ray electron and proton spectrum, and 
H.E.S.S. observes only 2900 6=*= with ( > 0.9 compared 
to 5020 protons. 

4. RESULTS 

A remarkable aspect of ACT searches for e^ anisotropy 
is that no dedicated observations are required in addition 
to archival data. As long as the target pulsars fall in re- 
gions of the sky such that, on average, the ACT spends 
a significant fraction of observing time in hemispheres 
containing and not containing the pulsars, a sizable frac- 
tion of the total archival observation time can be used 
to search for e anisotropics. The two pulsars we con- 
sider here (Monogem is at (l,b) = (201.1, -1-8.3), and 
Geminga is at (l,b) — (195.1, -1-4.3)) fall, for example, 
in a favorable part of the sky for H.E.S.S. observationfR 
In practice, the available observation time is the smalT 
est of the collective time the telescope spent observing in 
hemispheres containing and non-containing the relevant 
candidate pulsar. 

In order to predict the H.E.S.S. constraints which could 
be set using current data, we need to make a multitude 
of assumptions. We caution the Reader that several rele- 
vant figures have been approximated quite roughly, as an 
accurate calculation is impossible without further infor- 
mation about the H.E.S.S. telescope as well as its point- 
ing history. We first approximate the total available 
livetime of the instrument to be approximately 3000h 
for searches off of the Galactic plane (where the 7-ray 
background is much more significant), noting that the 
recently published H.E.S.S. line search included 1153h of 
livetime taken du ring studies of extragalac tic objects be- 
tween 2004-2007 ( [Abramowski et al.|2013| , and that ad- 
ditional surveys (such as those ot dwart spheroidal galax- 
ies) will also provide excellent targets for e* anisotropy 
searches. In order to calculate the effective area of the 



H.E.S.S. instrument, we take the model of Benbow ( 2005 
Fig 2a.), assuming an average zenith angle ot 45" , and we 
take the H.E.S.S. field of view to cover 3.8x10-^ sr. We 
additionally consider a future global observation time of 
5000h, including new data taken by the H.E.S.S. -II tele- 
scope. 

We note that this combination provides an effecti ve ex- 
posure at 340 GeV which exceeds that reported by |Aha-| 
ronian et al. (2009) by approximately a factor of 8, pre- 
sumably due to additional cuts regarding the nature of 
the electron shower and the removal of 7-ray point source 
contamination. We note that this mismatch is slightly 
wo rse when comparing th e 1 TeV calculation with that 



of (Aharonian et al. 20081, as additional cuts are made 
in this study. In what follows we degrade our calculated 
effective exposure by an ad hoc factor of 5. We note that 
this may be overly conservative, as cuts optimized for 
this study could allow for a much larger effective area, 
since this analysis is by its nature less susceptible to sys- 
tematic errors regarding proton contamination. 

In order to calculate the effect of the ( parameter on 
the separation of electrons and protons, we make a cut 

* We note that both pulsars are located relatively close to the 
Crab nebula, which is closely monitored by all ACTs. 
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Fig. 2. — The limits on the cosmic-ray e^ anisotropy one year of 
Fermi-LAT data (orange triangles) and those recently reported by 
AMS (cyan), as well as the predicted limits from 5 and 10 years of 
Fermi-LAT observations (orange solid and orange dashed), along 
with the predicted hmits from 3000 and 5000 hr of H.E.S.S. ob- 
servations (maroon solid and maroon dashed), as well as predicted 
limits from 1000 and 3000 hours of CTA observations (blue solid 
and blue dashed). These limits are compared with the predicted 
fluxes for models of the Geminga (black solid) and Monogem (red 
solid) pulsars which correctly explain the positron excess observed 
by AMS-02. We note that limits from the Fermi-LAT are techni- 
cally set based on a minimum energy E, rather than a traditional 
E dN/dE, a difference which is less important given the steeply 
falling e^ flux. 



of C > 0.9 J and use the values giv en by (Aharonian et al. 
2008 1 and Aharonian et al. ( 2009 ) to calculate the loss ot 
effective area after this cut is applied at energies of 1 TeV 
and 340 GeV respectively, which correspond to the low 
end of the energy range employed in each analysis. Thus, 
we assume a relative electron acceptance of 0.38 at en- 
ergies below 340 GeV and 0.77 at energies above 1 TeV, 
linearly interpolating between these values at intermedi- 
ate energies. Finally, using the observed proton back- 
ground at C > 0.9 for each observation, we assume an 
irreducible proton background which is 1.73 (0.95) times 
larger than the total e flux at energies of 340 GeV (1 
TeV), and again linearly interpolate between these val- 
ues. 

In order to compare our results with the projected lim- 
its from five years and ten years of Fermi-LAT dat a, we 
take the limits from one year of data obtained by |Ack- 
ermann et al. (2010b) and calculate a best fit power law 
to the 95% upper hmits of A < 0.30 (E / 1 TeVY'^^ 
and then reduce the limits by the square-root of the 
additional exposure time, in order to extrapolate opti- 
mistic results where the exclusion limits are dominated 
by statistics only. We note that this result may slightly 
improve, in light of new data-taking algorithms such as 
Pass-8, which are hkely to increase the effective area to 
electron showers (Atwood et al. 2013). We addition- 
ally compare our result with the 95%i confidence linrit 
of A < 0.036 measured by AMS-02, though we note that 
no energy scaling was provided with this value. 

In addition to these H.E.S.S. and Fermi-LAT ob- 
servations, we predict observations from the upcoming 
Cherenkov Telescope Array (CTA). While the parame- 



ters of CTA are unknown, for the purposes of repeata- 
bility we assume the folfowing educated guesses for the 
relevant parameters. We take the effective area to be 
an order of magnitude greater than the H.E.S.S. tele- 
scope, and the field of view to be a factor of (3/2)^ larger, 
corre s pondin g to the models suggested iu jCTA Consor-j 
Additionally, we note that the multiple 



tium 



^_^ ( |2011| ). 

telescopes of CTA greatly enhance its hadronic rejection 
capabilities, thus we add an ad hoc factor of two decrease 
in the relative hadronic flux for CTA observations. All 
other parameters (including the C cut and the factor of 
5 degradation in the effective exposure seen in H.E.S.S. 
observations) are left the same. 

In Figure [2] we show the current limits given by 1 year 
of Fermi-LAT data, the recently released AMS-02 lim- 
its on the e* anisotropy, as well as the projected lim- mum 7-ray flux of 1.33 x 10~^ TeV~^m~^s~^sr~^, which 



7 = 2.41 ( |Abdo et aL]|2010|. This implies that at en- 
ergies near 1 i'eV, the FGKB is subdominant to the e* 
contribution by a factor of ^-^250. Thus, we do not con- 
sider the effects of the EGRB on our results. 

A larger concern stems from Galactic diffuse 7-ray 
emission, which (due to our position in the Galaxy) has 
a dipole anisotropy of order unity. In order to ascertain 
the influence of this component on our results, we employ 
Galprop models of the Galactic diffuse emission, follow- 
ing the same parameters employed previously in the text. 
At high energies, we find the dominant source of diffuse 
7-ray emission to stem from tt" decays produced when 
energetic protons collide with Galactic gas. Eliminating 
observations within 5° of the Galactic plane, as well as a 
10° X 10° box around the Galactic center, we find a maxi- 



its from 3000h and 5000h of H.E.S.S. observations and 
from lOOOh and 3000h of CTA observations, compared 
to the projected anisotropics of the Geminga and Mono- 
gem pulsars with the same setup as described in Section 
[2] and shown in Fig. [l] Note that the observation times 
are meant here as total duration of the two ensembles 
of observations in the hemispheres towards and opposite 
the direction of the pulsar of interest. For example, for 
3000 hr, in the notation introduced in the previous sec- 
tion, the global duration of the ensemble of observations 
yielding Nf^^ is 1500 hr each. For all models we assume a 
total e flux given by th e best fitting power-law given by 



Aharonian et al.| ( |2009[ ). We find that current H.E.S.S. 
archival observations have the potential of observing the 
anisotropy induced by Monogem, while CTA observa- 
tions will be necessary in order to observe any anisotropy 
from Geminga. The predicted level of anisotropy from ei- 
ther pulsar under consideration here is fully compatible 
with the hmits from Fermi-LKAT and from AMS-02. 

5. THE DIFFUSE GAMMA-RAY BACKGROUND 

A serious concern in the identification of the e* 
anisotropy stems from the misidentification of the 7-ray 
background. Unlike hadronic showers, 7-ray and electron 
showers are nearly identical in nature, with the only ob- 
servable difference pertaining to a slightly different value 
of the reconstructed shower maximum (J^max)- How- 
ever, this parameter is very difficult to accurately mea- 
sure with ACTs, and thus only approximately 20% of the 
7-ray signal can be effectively eliminated before signifi- 
cantly cutt ing into the electron effective area ( Aharonian] 
et al.J[2008|). 



'I'here are two important sources of 7-ray contamina- 
tion that we need to consider here. The first corresponds 
to the extragalactic 7-ray background (EGRB), which 
is the the brightest component of the very high energy 
7-ray sky, and is thus the major contributor to the un- 
certainty in the e* spectrum. However, the EGRB is 
extremely isotropic, especially on large scales, and does 
not contribute to the e^ dipole anisotropy. The EGRB 
can thus be treated as a statistical background, as in 
the case of the hadronic contribution. In order to calcu- 
late the size of this effect, we extrapolate the calculated 
EGRB from Ferm-LAT observations, which indicate the 
EGRB to follow a power law 

dN/dE = k{E/lTeV)-'^, (13) 

with k = 3.3 X 10""^ TeV^^m-^s-isr-i and 



falls below the total lepton flux by a factor of 60. Ad- 
ditionally, by setting somewhat stronger constraints to 
limit the regions where Galactic gas may produce bright 
emission (setting |6| > 10° and eliminating a box of 20° 
around the Galactic center), we find a maximum 7-ray 
flux of 9.0 X 10"^ TeV~-^m~^s~-^sr~^, which is subdom- 
inant by more than two orders of magnitude. 

Since the diffuse 7-ray flux contributes with order 1 
anisotropics, its contribution to the overall anisotropy 
can still exceed the predicted e* anisotropics, which are 
likely to be less than 1%. However, we note that this sys- 
tematic error can be controlled in two independent ways. 
First, the dipole anisotropy stemming from 7-ray emis- 
sion will produce a peak in emission pointing towards the 
Galactic center, while the Monogem and Geminga pul- 
sars are both located in the Galactic anti-center region. 
While the pulsars may have moved consider ably since 



their e* flux was produced (Profumo 20121, it is un 



likely that either pulsar was located in a position aligned 
with the Galactic center. Secondly, we note that the 
morphology of tt*^ decay is relatively energy independent, 
and is extremely well mapped at energies below 100 GeV 
by the Fermi-LAT, thus a statistical comparison of any 
anisotropy with the Fermi-LAT diffuse 7-ray sky should 
immediately indicate, or rule out, a 7-ray origin. This 
method can additionally be used to remove portions of 
the e^ anisotropy which correlate strongly with locations 
of high diffuse 7-ray luminosity. 

6. SYSTEMATIC UNCERTAINTIES 

There are several systematic issues which further com- 
plicate the analysis presented here, and which will need 
to be carefully considered in any experimental effort to 
detect an electron anisotropy. First, unlike the Fermi- 
LAT, which achieves nearly uniform sky coverage, ACTs 
observe non-uniform patches of the sky. While no direct 
observations of candidate pulsars are necessary in order 
to perform the calculation suggested here, the ability for 
ACTs to determine the e^ anisotropy depends sensitively 
on the effective exposure of observations in the direction 
of each pulsar. For Geminga and Monogem, we note that 
their relative proximity (within 10° on the sky), and their 
location near the galactic anticenter, makes this obser- 
vation reasonable. 

Additionally, we note that our assumption of proton 
isotropy does not hold for the very precise measurements 
ACTs are capable of making. While the exact magni- 
tude of this anisotropy is difficult to ascertain, as mea- 



suremen ts by MILAGRO (|Abdoet aIJ[2009|) and ICE 



CUBE (Abbasi et al. 



2009[ |2UlU| )^neasured the pro 



ton anisotropy tor median proton energies of 6 TeV and 
20 TeV respectively. However, both experiments observe 
anisotropics on the order of 0.1%, of similar magnitude 
to the predicted lepton anisotropy. While a separation 
of these possible anisotropics requires careful experimen- 
tal consideration, some approaches exist which may help 
separate these components. For example, we note that 
the proton measurement should be exquisitely measur- 
able by ACTs, with almost no electron background con- 
tamination, simply by evaluating showers in ranges of 
C where the proton flux dominates by orders of mag- 
nitude. The resulting proton anisotropy can then be 
subtracted from the anisotropy of the electron measure- 
ment by employing a template fit. We note that the 
pulsar sources which we expect to dominate the electron 
anisotropy emit a negligible proton population, making 
the measured anisotropics largely uncorrelated. We also 
note that the proton anisotropy is observed to exist in lo- 
calized hotspots, which may be separable from the dipole 
dominated electron anisotropy our models suggest. 

Finally, we note that anisotropics can be induced due 
to instrumental systematics, such as uncertainties in the 
effective acceptance could artificially induce an instru- 
mental sensitivity. We note that while overall errors in 
the calculated effective area are immaterial to calcula- 
tions of the anisotropy, both optical degradations (espe- 
cially in the photo-multiplier tubes) and night-to-night 
optical variations (based on the the background starlight 
and ambient temperature) can have a large effect on 
the measured anisotropy. Again, controlling these sys- 
tematics requires careful experimental work in evaluat- 
ing the observed electron spectrum under variations of 
many known parameters, and we can only suggest mech- 
anisms which may make this measurement possible. We 
note that the event-shuffling technique employed by the 
Fermi-LAT collaboration to study anisotropics can, in 
principle also be employed with AC T observations in 
order to test for systematic biases (Ackermann et al.| 
|2010b). Furthermore, the large statistics m the proton 
anisotropy detected above can be used to independently 
test the systematic uncertainties in the instrumental ef- 
fective area. 

7. DISCUSSION AND CONCLUSIONS 

We showed that the recent measurement of the 
positron fraction by AMS-02 is in excellent agreement 
with the simplest models of a diffuse cosmic-ray back- 
ground with contributions from a single, nearby pulsar 
emitting e pairs with a power law spectrum with spec- 
tral index close to 2. Since this stands as the minimal 
explanation which fits all current data, arguments along 



Occam's Razor would make it a preferable model. Sec- 
ondly, we showed that this scenario is differentiable from 
both models of dark matter annihilation, as well as from 
models invoking a multitude of pulsars, based on the in- 
duced anisotropy in the cosmic-ray e* arrival directions 
that would stem from one single (or very few) bright, 
nearby source. Finally, we have proposed a novel method 
to observe the predicted e* anisotropy, using the wide 
collecting area of ACTs and large ensembles of archival 
observations in order to collect a substantial population 
of cosmic-ray e^ from all directions in the sky, and ar- 
gued that this method is capable of setting the strongest 
constraints on the cosmic-ray e^ anisotropy. 

We note that, under certain assumptions on the dark 
matter mass and annihilation final states, relatively large 
anisotropics can be produced even in the context of 
dark matter annihilation as t he source o f the anomalous 
positron fraction (BorricUo et al. 2010| ). However, we 
note two possible complications for dark matter models 
contain ing significant anisotro pics. First, the predictions 



of Ref. Borriello et al. (20101 are significantly affected 



by AMS-02 results, since very special dark matter final 
states and masses are comp atible with the high-stati stics 
AMS-02 data (see e.g. Ref. |Cholis fc Hooper| ( [20131 )); it 
is unclear whether th e large levels of anisotropy found in 
Borriello et"al] ( |2010[ ) are still viable after the AMS-02 re- 
sults. Second, the largest anisotropies from dark matter 
stem from a highly clum py dark matter densi ty distribu- 
tion, as discussed in Ref. Borriello et al.| (2010). However, 
a bright local clump would h ave been ob ; 
rays, as pointed out e.g. Ref.|Brunet al 



s erved in gamma 
(|2009|); it is thus 



unclear whether any large anisotropy would be actually 
compatible with gamma-ray observations. 

Finally, while this work has focused on prospects for 
H.E.S.S., the position of VERITAS in the northern hemi- 
sphere makes it an equally good telescope with which to 
measure the e^ anisotropy independently of the H.E.S.S. 
results. While systematic effects may greatly degrade the 
sensitivity of each instrument towards the detection of 
electron anisotropies, we note that ACTs arc predicted 
to set constraints on the TeV scale electron anisotropy 
which are nearly two orders of magnitude stronger than 
the predicted 10-year constraints from the Fermi-LAT. 
This means that even if systematic issues greatly degrade 
the limits which can be set by ACT instruments, they are 
still likely to probe interesting regions of the parameter 
space. 

We would like to thank Amy Furniss for valuable 
discussions on the capabilities of ACT telescopes. This 
work is partly supported by the Department of Energy 
under contract DE-FG02-04ER41286 
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